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Mathematical modelling of metabolism 
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Mathematical models gf the cellular metabolism have a special 
interest within biotechnology. Many different kinds of 
commercially important products are derived from the cell 
factory, and metabolic engineering can be applied to improve 
existing production processes, as well es to make new 
processes available, Both stoichiometric and kinetic model* 
have been used to investigate the metabolism, which has 
resulted in defining the optimal fermentation conditions, ae well 
as in directing the genetic changes to be introduced in order 
to obtain a good producer strain or cell line. With the 
increasing availability of genomic Information and powerful 
analytical techniques, mathematical models also serve as a tool 
for understanding the cellular metabolism and physiology. 
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Abbreviation 

M FA metabolic flux analysis 

Introduction 

Mathematical modelling is a very powerful tool in physics, 
chemistry, and engineering for interpretation and predic- 
tion of natural phenomena and experimental results [1]. A 
mathematical model is always a simplification of the actu- 
al phenomenon and it is therefore possible to establish 
different mathematical models for che same phenomenon, 
depending on the objectives of the model and the avail- 
able measurements (2). 

With the explosion in experimental data within biology, 
there arc many attempts to develop mathematical models for 
description of cellular functions — either overall function or 
the function of individual cellular processes. Thus, in the 
field of functional genomics, bioinformatic* play* a promi- 
nent role in assigning function to orphan genes. Within 
biotechnology there is especially a focus on the . cellular 
metabolism, as this may be exploited for the production of 
compounds that might find application as materials, pharma- 
ceuticals, food additives, and so on. Besides featuring: very 
complex networks, with interconnecting pathways thac con- 
sist of hundreds of reactions, the metabolism of a living cell 
is also subject to control and regulatory mechanisms. These 
regulatory mechanisms are not completely elucidated and 
arc therefore very difficult to quantify. Thus, the establish- 
ment of fully mechanistic models to describe the cellular 
behaviour in terms of its metabolism is, at least to date, not 



possible and all models are therefore based on simplifica- 
tions. In this context, mathematical models play a crucial role 
in hypotheses testing, that is, they can^ serve as a guide to 
choosing among different possible regulatory structures for a 
specific cellular process. 

In traditional studies of fermentation processes, extracel- 
lular metabolites (such as substrates and products) have 
been measured, as well as the biomass concentration. 
The models that can be formulated based on these type 
of measurements are, however, highly unstructured and 
their application to the interpretation! of cellular physiol- 
ogy and prediction of cellular behaviour under different 
cultivation conditions i* quite Hmhjcd. Following the 
advances in analytical techniques, measurement of the 
concentrations of intracellular metabolites and the activ- 
ities of intracellular enzymes ha* allowed for the 
formulation of more structured models, that have 
increased the possibility for interpretation and predic- 
tion of cell physiology. Besides this, the increase in 
speed and power of computers has made it possible to 
solve the highly non-linear models that arose from the 
inclusion of more structure in the models. More recent- 
ly, the development of new powerful analytical 
techniques such as DNA arrays, 2D-p\ electrophoresis, 
and mass spectrometry has enabled a very detailed 
analysis of cellular function that might form the basis for 
a new category of mathematical models that describe the 
overall cellular function. 

In this review, we focus on mathemadcal;mode!s that describe 
the cellular metabolism, as these models; play a central role in 
the rapid developing field of metabolic ^engineering [3,4], In 
our discussion, we group the models according to their struc- 
ture; firstly, stoichiometric models, which are based on the 
time invariant characteristics of metabolic networks; and sec- 
ondly, kinetic models, which are usually based on both 
stoichiometry and enzyme or microbial kinetics. 

Stoichiometric models 

Metabolic flux analysis (MFA) has been widely used for the 
quantification of the intracellular fluxes in the central metab- 
olism of bacterial, yeast, filamentous fungi and animal cells 
(Figure 1). In MFA, mass balances ovct all the intracellular, 
metabolites are used to calculate the fluxes through the dif- 
ferent branches of the network. Hereby it is possible to get a 
snapshot of the metabolism under a particular condition. 
The fluxes can be calculated by combining measurements of 
a few fluxes either with linear algebra Or linear optimisation 
[5,6]. More recently, the use of labelled substrates combined 
with the measurement of the labelling! state of intracellular 
metabolites, either by NMR [7,8*] or by gas chromatogra- 
phy/mass spectrometry (OC-MS) [9Ji has been used to 
estimate the fluxes. When the material balances used in 
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Principles o1 stoichiometric modelling. Firstly, 
a stoichiometric matrix is defined that 
appropriately describes the metabolism under 
investigation, baaed on available biochemical 
information. Secondly, the stoichiometric 
matrix is multiplied by the ao-called vector of 
reaction rates, defining the mass balance 
equation. Finally, algebraic manipulation Of the 
stoichiometric matrix can be carriod out by 
different means, depending on the objective 
of tho analysis, (a) It is possible to impose 
some constraints (e.g. by measuring some 
fluxes) SO that the system becomes 
determined and can be solved by simple linear 
algebra. In this case, the analysis is usually 
called flux balancing. (b> If further constraints 
are imposed by measuring the labelling stato 
of some metabolites, however, it may be 
possible to combine flux balancing with 
labelling balancing and the system has to be 
served numerically. In this case, a better 
estimation of the fluxes is attained, (c) If the 
system cannot be constrained to a 
deter mined one, linear optimisation can be 
applied to find the maximum Of minimum of a 
suitable objective function, (d) Instead of 
calculating particular solutions, it is possible 
to gain aome insight into pathway structures 
within a given metabolic network by applying 
convex analysis to determine the so -called 
elementary flux modeB or (e) by calculating 
biochemically meaningful base vectors for the 
null space of tho stoichiometric matrix. F6P, 
fructose-6-phoSphate: G6P, glucoae-B- 
phosphate; GLC, glucose. 
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craditional MFA arc combined with balances of the labelling 
pattern of the metabolites, chc models become non-linear 
{10-13 J 41- The additional information s applied by mea- 
surements of the labelling pattern of the metabolites do, 
however, allow for a more reliable estimation of the fluxes, as 
well as the analysis of the pathway scructurcs and possible 
reversibilities. These new features make it possible co spec- 
ulate about some key points in cell metabolism, such as 
pathway identification and compartmcntation of enzymes 
and metabolites, which lead to the coining of the ccrm 'meta- 
bolic network analysis' (MNA) [14"J. 

Besides being applied to the investigation of cells under dif- 
ferent environmental conditions and for studying 
different mutants of a particular cell (see Table 1 for some 
recent MFA applications), these stoichiometric models have 



also been applied with the aim of predicting the 
gcnotypc-phcnotypc relationship, in order to fill in the gap 
that exists between DNA sequence data^and functional infor- 
mation. Schilling etdl. [15'J have constructed a stoichiometric 
model describing ail known reactions in EsrAtrirJiut coR^ and 
have analysed the phenotypes of different deletion mutants, 
using linear optimisation. The predictions made by using this 
approach were in agreement with 60 of the 66 mutants exam- 
ined. Furthermore, it was proposed that by identifying a 
proper set of vectors that span the null space of the stoichio- 
metric matrix or by determining the so-called elementary flux 
modes by convex analysis, it is possible to assess all the capa- 
bilities of a metabolic genotype (16,17**j (Figure 1). 

Stoichiometric models arc clearly very powerful, but 
their main drawback is the limited predictive power, 
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Tablet 



Soma recent applications of stoldilo metric models. 



Syslem investigated 



Kind of approach Kind of analysis 



Main achievements 



Serine alkaline proteaaa 
production by B. lichenihrmia 

£. co ft mutants lacking the 
acetate production pathway 
(aefc/4-pffl genea) 



£ coA' expressing tho 
acatolactate synthase 
gBne of ft subtitis 

£ coli expressing the phb 
oporon Of ft eut/vpha 



Pyruvate production by 
T. glabmta 



Production of IgG against 
human fibronectin by murine 
hybridoma cells 

Wild-type and recombinant 
Z. moo/7/s fsrmantation on 
glucose, fructose and xyloaa 



MFA + linear 
optimlaation 

MFA 



MFA + linear 
optimisation 



MFA 



Refercnco 



MFA 



MFA + 
labelling experiments 



Production of PHB in mlxod 
cultures of L dclbrvtckti and 
A. eutrophua 



Citrate-glucose cometabolism 
In a tubtilis 



Enhanced biological phosphorus 
removal (EBPR) 



MFA 



MFA 



MFA + 
linear 
optimisation 



(a) Design of an oxygen transfer strategy (42) 
in the bioreactor for SAP production 

(not verified oxperimentaOy).. 

(b) Mutations in both vckA-fit* ! [43] 
and nuo genea (encodes for a 

NAOHrubtquinone dehydrogenase) 
are required for reduction of thojflux 
through the pyruvalo-formate-lyase 
reaction. 

(b) Reduction of acetate formation [44) 

and increase of acetoin formation, 
which is less toxic to growth and 
heterologous protein eyntheelk 

(b) In order to maximise the yield (45] 

of PHB, tho availability both of acetyl 
CoA and NAOPH have to be 
increased (experimental reeuljs 
were not in agreement with; 
MFA predictions). 

(a) Pyruvato production ia mainty [46] 

affected by the fluxes through pyruvate 
dehydrogenase and pyruvate 
da csrboxy lase reactions. A thiamine 
feeding strategy and the influence of 

dissolved oxygen on pyruvate 
accumulation were also invesligiled. 

Interpretation of steady-state [47] 
multiplicity (different cell concentrations 
at the same dilution rato) as a different 
utilisation of pyruvato in the TCA Cycle. 

Synthesis of ribo*e-5-phoephate Occurs [48] 
malhty via the reversible transJcet^lase 
reaction. Evidence of reversibjo 
operation of riboso-5-phosptane 
isomerase, phoephoglucoae isomtarasa 
and ribulose-5-phoephate epimeraae. 
Identification of the xyMokinaae ehsyme 

activity as the rate controlling step 
for ethanol production in recombinant 
Z mobifis. 

(a) NADPH generated in the isocitiate [49] 

dehydrogenase catalysed reaction in 
A eutrophua is mainly used for ot-ketygfutarate 
conversion into glutamale when 
NH 3 is abundant, whereas K is usjod for 
PHB formation whon NH a j 
concentration decreases. ■ 

(a) ' Cofeeding of citrate with glucose 1 [50] 

reduces acid formation and increases 
cell yield on carbon source due io an 
attenuation of pyruvate kinase flux. 

(a) Probably the first attempt to appr> MFA [51] 

to an undefined, mired culture. 
Degradation of polymers and uptake 
or release of phosphate and acetate in an 

EBPR system could be predicted by 
the model; however, the results differed 
up to 55% from experimental data. 



(a) 



(a) 
(b) 



(a) Same cells under different environmental conditions, (b) Different mutants under the same environmental conditions. 
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which is due to the lack of regulatory information in the 
model formulation. 

Kinetic models 

When detailed information is available about the kinetics 
of specific cellular processes (e.g. enzyme-catalysed reac- 
tions, protein-protein interactions, or protcin-DNA 
binding) it is possible to describe the dynamics of these 
processes by combining kinetics with the known stoi- 
chiomctry of metabolic pathways (Figure 2). Rizzi et aL 
[18] have applied this methodology to model glycolysis in 
Saccharomyces ccrevixiat- The model includes reactions of 
the Embden-Meyerhof-Parnas (EMP) pathway, tricar- 
boxylic acid (TCA) cycle, glyoxylate cycle, and respiratory 
chain (22 material balances around metabolites and 23 rate 
equations for enzymes have been included) and Has been 
successfully applied to predict the levels of intracellular 
and extracellular metabolites after a glucose puUe in a con- 
tinuous culture of S. cercvisiaz. The model has only been 
applied to a 120 s ti mescal e, however, and phenomena 
such as enzyme synthesis and degradation — which would 
be important in larger timescales — have not been includ- 
ed. Besides this, regulation has only been included at the 
individual enzyme level and all the kinetic parameters for 
enzymatic activities arc taken from in vitro experiments, 
which probably do not reflect the in vivo situation. Vascghi 
et al, [19*') have applied the same strategy to investigate 
the pentose phosphate (PP) pathway in S, cenvisiae. The 
concentrations of some intracellular metabolites, such as 
glucose-6-phosphatc and 6-phosphoglueonate, as well as of 
the co-enxyme NADPH predicted by the model were in 
agreement with experimental results. By using the model, 
some kinetic rate equations and parameters were identi- 
fied — as they could not be found directly in the 
literature — and it was suggested that the split ratio 
between the PP and the EMP pathways is related to the 
intracellular concentration of MgATP*-. 

Kinetic modelling has also been applied to investigation of 
the penicillin biosynthctic pathway in Praia 'ilium chryso- 
gfmtrn [20-22]. The model was used for calculating the 
fluxes through this pathway (enzyme kinetics for 10 reac- 
tions were included), as well as the concentrations of the 
involved metabolites, which were in agreement with 
experimental results. In these studies [20-22], the models 
have been applied to analysis of the flux control in the 
pathway, which demonstrates the use of kinetic models in 
the field of metabolic engineering. Van Ricl etaL [23J have 
established a kinetic model for the central nitrogen metab- 
olism in S. terevisiat^ including variables that account for 
regulatory aspects. These variables represent the concen- 
tration of so-called regulators and may indicate a 
way — though non-mccanistic — of accounting for the 
control aspects of metabolism. 

Another way of accounting for the regulatory aspects of 
cell metabolism is by applying cybernetic principles 
[24,25]. In this case, cybernetic variables are introduced 
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Principlea of kinetic modelling. In this Case, a simple reaction network 
etoichiornetry (a) is UBed in combination with simple Michaete-Menten 
type enzyme kinelica (b) to yield the mass balance over metabolite M 2 , 
as shown (c>. If the parameters and all tine Michaelia constant 
for reaction r^) are known, ae well aa the concentration of metabolite 
W\ y (which may be a substrate), it is possible to calculate the timo 
profile of metabolite M 2 . 



into a kinetic model with the aim of substituting the 
unknown mechanistic details of the cell regulatory archi- 
tecture by an objective function by Supposing that the 
metabolism of a cell operates with a specific overall goal, 
such as the optimisation of growth. This approach has 
been used to predict the increase in jflux towards threo- 
nine formation in Coryncbacterium laetofcnnetitum^ as a 
consequence of genetic modifications! [26]. A further way 
of including regulatory aspects in a mathematical model is 
the approach presented by Hatzimanikatis tt al. [27], 
which formulates the model as a 'mixed integer linear pro- 
gramming 1 (MILP) optimisation problfem. In this case, the 
presence or absence of different possible regulatory loops 
are represented by discrete variables/ This approach has 
been applied to the investigation of genetic alterations in 
recombinant high-cthanol-producing E. cali [28*]. By 
using a (log)linear approximation for a non-linear model, 
the authors verified a qualitative agreement between the 
predictions of the model and experimental results. When 
mechanistic details arc missing, ncuraf networks [29,30] or 
fu/^y logic-based models [31,32] can Be used with the aim 
of simulating metabolic behaviour; however, the amount 
of experimental data required to generate these models is 
usually very high and also less insightiinto the underlying 
mechanisms is gained. To illustrate the difficulty in mech- 
anistically modelling particular aspects of metabolism, 
Wong et aL [33] presented a model for the lac opcron — 
including aspects such as catabolitc repression and induc- 
er exclusion — that required several assumptions to be 
made. Thus, even in this very well studied system there 
are still some mechanistic aspects that have nnt been elu- 
cidated. Some other recent applications of kinetic models 
are listed in Table 2. 
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Table 2 



Some ntcent applications of kinetic model*. 


System investigated 


rstna or approacn 


Main ar.hiflvflmftnta 


Reference 


PHB synthesis pathway in 
A outrophas 


Investigation of the sol b 
PHB pathway by using 
onzyme kinetics 


It is shown that the use of oversimplified 
kinetic expressions, e.g. MicheWis- 
Mcnten type, is inappropriate for describing 
PH 3 synthesis, and also that parameters 
taken from literature should have boon 
determined under the appropriate 
physiological conditions. Simulations give 
support to AcCoAiCoA and NAQHPK: 
NADP+ ratios exerting regulation on 
thiclase and reductase catatyeed 
reactions, respectively. 


152] 


Biological phosphorous 
removal by Acinetobacter so. 


Investigation of batch 
cultivations by structuring 
bbmasa into three compartments 


Tho model was successfully i&sed 

lu pro QIC I lilO IRuISWIH FWpLifipf* u ' 

Aa'netobacter sp. cultures from fcerobic 
to anaerobic environments, and vice versa. 


(53] 


Hyaluronic acid fermentation by 
S. zoocpidcmtCU* 


Investigation of batch and 
fod-batch cultures by structuring 
biomass into two compartments 


ft was possible to predict the 
bio mass and hyaluronic acid iima 
profiles in aerobic and anaerobic 

culture* of 5, zocGptdemicba. 


CS43 


Production of cellutaae 
by r /-eese/ 


Investigation of batch cultivations 
by three types Df models: unstructured, 
Structured into three compartments, and 
using neural network parameter function. 


The prediction capacity and the process 

optimisation results using the neural- 
network -based model were higher when 

compared to the unstructured and 
structured models tested. However, both 
the neurar-network-baaed 
model and the unstructured niode] 
do not give any insight into the system's 

underlying mechanisms., 

_ } 


[55] 



If information on all this enzymatic reactions (both equa- 
tions and parameters) of the whole metabolism of an 
organism were available, ic would in principle be possible 
to apply detailed modelling to interpret experimental 
results and predict the dynamic behaviour of cells when 
subjected to a determined shift in the environmental con- 
ditions and the consequences of specific genetic changes. 
The limitation of this kind of approach, however, is the 
lack of information on general regulatory aspects, as well as 
the fact that normally only in virn? parameters arc available 
for the enzyme kinetics, in vivo perturbation experiments 
could be a way of fine-tuning these parameters, as illus- 
trated by Rim** a/. [18]. 

Mathematical modelling of regulatory phenomena, such as 
heat-shock regulation [34], synergistic cukaryotic gene 
activation [35 s ] v regulation of the Gl— >S transition of 
cukaryotic cells [36J, signal transduction pathways [37]» 
and gene expression in a complex rcgulon (38*), can serve 
as a tool — which should be used together with expert 
mental investigation — in clarifying the molecular 
mechanisms behind these phenomena. This may lead to 
combining modelling of the metabolism with modelling of 
signal transduction pathways. 

Wich the aim of trying to interpret die genotypc-pheno- 
type relationship, Hauimanikatis et a/. [39) propose the 



dynamic analysis of the protein levels of an organism 
instead of using a set of metabolic fluxes as the output of a 
non-dynamic model. As it is not possible to correlate 
mRNA and protein levels directly [4Q], the sole informa- 
tion on gene expression generated by hybridisation on a 
solid surface (such an from DNA mic roar rays) cannot 
account for the functional assign mcjnt of the genome 
sequence of an organism. In order to illustrate the need of 
combining mRNA and protein levels, iHatzimanikatis and 
Lcc [41 M ] applied a continuous dynamic model to 
describe arcadian rhythm icity that agrees with experimen- 
tal observations. In spite of being incipient as a result of 
limitations that still exist in protcomc analysis, this kind of 
functional modelling is very promising as it really focuses 
on the final product of the gene, which is the protein. 

Conclusions and future perspectives 

Modelling of metabolism is either predicting the cellular 
behaviour under a different cnvironmcbl conditinn, given a 
fixed genetic background, or predicting the cellular behav- 
iour under a different genetic background, given fixed 
environmental conditions. In both cases, the output is the 
phenotype, which may be described by metabolic fluxes 
and/or the concentrations of transcripts*! proteins or metabol- 
ic intermediates. Describing metabolism using a 
combination of these different characterisation forms — both 
at dynamic and steady-state conditions'— will be important 
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not only for optimising biotechnological processes, but also 
in the field of functional genomics. The exercise of mathe- 
matically modelling metabolism and the involved regulatory 
phenomena perscc** certainty serve as a tool In the elucida- 
tion of the cells' control architectures. The limitation of all 
approaches is still how to include regulatory and control 
aspects of metabolism, which are either difficult to describe 
or still lack biochemical characterisation. In light of this, 
dynamic models arc not so widespread in literature, and are 
usually applied to the analysis of relatively small parts of the 
metabolism, compared to stoichiometric models, which can 
describe the complete metabolic network. 
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